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ABSTRACT 

We present the discovery of optical/X-ray flux correlations on rapid timescales in the low/hard 
state of the Galactic black hole GX 339-4. The source had recently emerged from out- 
burst and was associated with a relatively-faint counterpart with mag Vk.\1 . The optical 
(VLr/ULTRACAM) and X-ray (RXTEIVCA) data show a clear positive cross-correlation 
function (CCF) signal, with the optical peak lagging X-rays by ^ 150 ms, preceded by a shal- 
low rise and followed by a steep decline along with broad anti-correlation dips. Examination 
of the light curves shows that the main CCF features are reproduced in superpositions of flares 
and dips. The CCF peak is narrow and the X-ray auto-correlation function (ACF) is broader 
than the optical ACF, arguing against reprocessing as the origin for the rapid optical emission. 
X-ray flaring is associated with spectral hardening, but no corresponding changes are detected 
around optical peaks and dips. The variability may be explained in the context of synchrotron 
emission with interaction between a jet and a corona. The complex CCF structure in GX 339^ 
has similarities to that of another remarkable X-ray binary XTE Jl 118-1-480, in spite of show- 
ing a weaker maximum strength. Such simultaneous multi-wavelength, rapid timing studies 
provide key constraints for modeling the inner regions of accreting stellar sources. 

Key words: accretion: stars - individual: GX339^ - stars: X-rays: binaries - stars: optical: 
variable - black holes 



1 INTRODUCTION 

GX 339-4 is an X-ray binary system hosting one of the most 
promising Galactic bla ck hole (BH) can didates, with a mass 
function of '^ 6 Mq JHvnes et al.l 120030 . It is classified as a 
micro-quasar, based on tight radio/X-ray correlations found over 
a larg e dynamic range in source luminosity (e.g. iGallo et al.l 
l2003h . E xtensive X-ray studies have revealed a rich timing struc- 
ture (e.g. 'Punn et al J20o3:lBelloni et alj2005l ; lNowak et al.l 19991 : 
[Miya moto et al. 199 1|). In its low state, the source is typically as- 
sociated with an optically-bright counterpart, during which flick- 
ering on timescales as short as 10 ms has been observed. Rare 
transitions to the X-ray-off state have enabled plac ing some con- 
straints on the much fainter optical companion star JShahbaz et al.l 



* Based on observations carried out in ESO programmes 079.D-0535 and 
279.D-5021, and RXTE Proposal number 93119. 



1 200 11) . And during X-ray outbursts, the optical and X-ray fluxes 
display complex correlations as well as anti-correlations over ex- 
tended timescales (e.g. Makishima et al. 1986; Russell et al. 200q). 

Rapid timing observations simultaneous over a broad energy 
range remain as one of the lesser-explored aspects of GX 339^. 
The last published multi-wavelength studies to probe sub-second 
timescales were carried out over 20 years ago (but see Spruit et 
al. in prep.). Besides high state observations by Makishima et al.L 
lMotchetal.1 ( Il98lll983lll985h studied the source during a low- to 
high-state transition, and found an anti-correlated cross-correlation 
function (CCF) signal with optical leading X-rays by a few sec- 
onds. Although there was an indication in the low-frequency light 
curves that the optical and 13-20 keV X-ray fluxes were correlated 
at a significance of about 98 per cent, no obvious high frequency 
(< 20 s) correlation was uncovered at any of the selected energies. 

The Galactic-halo X-ray binary XTE J1118-H480 has recently 
been a target of extensive, simultaneous multi-wavelength cov- 
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erage. Among its interesti ng properties, a complex optical/X-ray 
CCF has been found (e.g. iKanbach et ai]|200ll : ISpruit & KanbachI 
l2002h . including positive and negative correlation components and 
an optical-vs. -X-ray peak lag of ~ 0.5 s. Com plex CCFs are now 
being found in other sources as well jDurant et al. 2008 ). These 
provide important time domain constraints for physical models. 

In this Letter, we present the first simultaneous optical/X-ray 
timing analysis of GX 339^ on rapid timescales of ~ 50 — 130 ms 
in a low/hard state associated with a relatively-faint optical coun- 
terpart. A clear cross-correlation signal is detected with the optical 
peak lagging X-rays by ~ 150 ms. The CCF has a complex pattern 
with some similarities as well as differences to the CCF found for 
XTE J1118-I-480, and is likely to be result of distinct interacting ac- 
cretion/ejection components. Full details of the timing and spectral 
analysis are presented in a forthcoming paper. 



2 OBSERVATIONS 



The triple-beam optical camera ULTRACAM (Dhillon et al. 2007), 
capable of high-speed photometry at up to 500 Hz, was mounted on 
the Very Large Telescope (VLT) as a visitor instrument during Jun 
2007. We carried out three 1 h long observations of GX 339-4 on 
alternate nights of UT Jun 14th, 16th and 18th (hereafter referred to 
as Nights 3, 2 and 1 respectively in order of improving weather), si- 
multaneously with the Rossi X-ray Timing Explorer (RXTE) satel- 
lite. Only the final night [Night 1] was photometric, while the other 
two had variable transparency. Night 3 being the worst. This period 
fell a few weeks after the source had returned to the low/hard state 
following a large X-ray outburst (Kalemci et al. 2007) and also co- 
incided with rising optical emission CBuxton & Bailva.2007.) . Our 
choice of time resolution was governed by the need to obtain a 
good signahnoise under prevailing atmospheric conditions. The fi- 
nal values used were ~ 50, 133 and 136 ms on Nights 1, 2 and 
3 respectively. Data calibration and relative photometry (with re- 
spect to a brighter comparison star observed simultaneously) was 
carried out with the ULTRACAM pipeline v. 8.1.1. Three-filter si- 
multaneous observations are possible and we used u' , g' and r', 
but much longer exposures were required in u' to obtain compa- 
rable signahnoise; consequently, the u' data is not considered fur- 
ther in this Letter. Optical spectro-photometry carried out with the 



VLT/F0RS2 instrument on Night 1 gives F^°°°^ 



erg 



^ (Vvcga ~ 17). Correction for Galactic extinction 

of ylv~3.3 im plies XLf°°^ » 4.6 x 10^^(d/8 kpc)^ erg s,"^ 
(distance from Izdziarski et al.ll2004h . 

RXTE observed the target in its canonical GoodXenon and 
Standard PCA modes, and in a 32-s on-off rocking mode with 
HEXTE cluster 1. Recommended HEADAS v. 6.4 procedures were 
followed for data reduction and extraction of light-curves and spec- 
tra, including the latest calibration information and background 
model corrections. A simple hard power-law with photon-index 
r = 1.65 ± 0.02 provided a statistically acceptable fit to the spec- 
tra for energies of 3-200 keV. The source had a flux F2-10 ~ 
1.6 X 10"^° erg s'^ cm-^ => L2-10 = 1.2 x 10^^ erg s"\ im- 
plying an optical(V^):X-ray(2 — 10 keV) luminosity ratio of about 
40 per cent. The low flux and h ard power-law are cha racteristic of 
the source in the low/hard state. iTomsick et alj ( 120081) find that the 
X-ray flux reached a minimum during our observation period. 

The source showed a high fractional variability amplitude in 
X-rays, approaching 50 per cent in the full-band PCA energy range 
(above the value expected from Poisson fluctuations) on all nights. 



In the optical, the net rms variability over the full light curves was 
« 13 and 15 per cent in the g' and r' filters, respectively. 



3 RESULTS : THE CROSS-CORRELATION FUNCTION 

The net optical and X-ray light curves were translated to a common 
Barycentric frame, and cross-correlated on the fastest optical time 
resolutions available on each night. The absolute and relative tim- 
ing accuracies of ULTRACAM are '--^ 1 ms and 50 /is, respectively 
( Dhillon et al. 20 07). much better than the smallest timescales in 
our light curves. The main result of our work is shown in Fig.[T] The 
CCF shows a single, significant peak at an optical lag of ~ 150 ms. 
The peak itself has a narrow core, with a shallow rise from -^ — 1.5 
s to s, and a steep decline from 150 ms to ~ 0.5 s. Weaker, but 
significant anti-correlation troughs are centred at ~ —4 s and 1 s. 
Each of these structures is visible in all the observations (in spite 
of some clear inter-night variation), suggesting that each of them 
is real. The peak narrowness and position are constant between the 
three nights, within the '--^ 50 — 130 ms resolution available. The g' 
data shows a very similar CCF to the r'/X-ray one presented. 

The asymmetric shallow rise and steep declin e of the CCF is 
clearly reminiscent of that seen in XTE Jl 1 18-1-480 ( IKanbach et al.l 
2001), but 'mirror-imaged' about a vertical axis and shifted to a lag 
of 0.15 s. This lag was confirmed by constructing averaged optical 
light curves around micro- (local) flares and dips selected in the 
full-band X-ray data. Fig. |2] shows the resultant optical (as well as 
2-5 and 5-20 keV X-ray) superpositions around several hundred 
X-ray peaks and dips. An optical extremum appears at ~ 150 ms 
lag on all nights (though only the best weather Nights 1 and 2 are 
shown). Furthermore, the peaks light curve clearly shows lower- 
than-average troughs, as well as a preceding rise of the optical, all 
matching the CCF within t=±2 s. There are also some matches be- 
yond this range, including the peaks local maximum at t=+10 s. 
But strong intrinsic (not Poisson) variability dominates, especially 
in the higher-resolution Night 1 data. This suggests that more flar- 
ing is present on timescales smaller than those probed by us. 

There may be sev eral reasons why the p r evious observations 
of iMotch ( 



•e may be sev eral reasons why the p r evious observations 
letal] h983l) and iMakishima eTal] jl98d) did not find 
a positive CCF signal. Firstly, the CCF strength itself probably 
evolves between the different states probed: a high-state during the 
observation of Makishima et al. and an optically-bright [V ~ 15.4] 
low state in iMotch et al J During our observations, the CCF in Fig.[T] 
has a peak strength of only ~ 0. 1; a perfect correlation would show 
a peak of 1 (for comparison, XTE Jl 1 18-1-480 has a peak of 0.4). 
It is also possible that the length of the simult aneous optical/X- 
ray observation (^ 100 s) available to lMotch et al. ( 1983i,) was too 
short to reveal any positive CCF components present. 



4 DISCUSSION 

The main question that we wish to address is the origin 
of the rapidly variable optical power and the complex cross- 
correlation. Published models for the XTE Jl 118-1-480 CCF 
include, among others: iMerloni et alj 1 20001), w ho invoke 
magnetically-dominated core 



Esin et al.l ( 1200 ll . 



a 
dominant 



advection-dominated flow [AD AF] with additional synchrotron), 
iMarkoff et al] ( 1200 iL a pure iet). lMalzac et alj(l2004 feedback in 
a common jet-l-corona reservoir) and I Yuan et alj ( 120051 an ADAF 
and a jet dominating at different energies). Most can explain the 
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Figure 1. The r' vs. X-ray full-band PCA cross-correlation for data on all three nights, obtained from simultaneous light curve sections of 60 s length. A 
positive delay (in this case peaked at Ri 150 ms) implies that optical lags X-rays. The inset shows a zoom-in average CCF of Nights 1 and 2 interpolated onto 
the fastest timescale of 50 ms in order to clearly illustrate the delay. The shaded region is the average scatter computed in an ensemble of light curve sections. 




Figure 2. Averaged 2-5 keV, 5-20 keV and optical (r') light curves around full-band PCA X-ray flares (left) and dips (right). Optical hght curves from Nights 
1 and 2 are shown in r ed and blue; averages from both nights are shown in black. For making this plot, full-band extrema are selected according to the method 
of lMalzac et al.l 120031 cf. their Fig. 9). A flare (or dip) must be at least f=2 times above (or 1// times below) the local X-ray mean in a running tm=32 s long 
section, and is also required to be the local extremum within a contiguous segment of ±tp=8 s (this effectively selects significant flares only). Corresponding 
light curves sections in other bands are then normalized to their local means before being averaged. EiTor bars show typical Poisson uncertainties. 
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Figure 3. The X-ray (black) and optical (red dashed) auto con'elation func- 
tions (ACF) computed from the highest time resolution (50 ms) light curves 
on Night 1 . The X-ray ACF is for the full-band PCA, and the optical refers 
to the r' filter. Both are corrected for Poisson noise (dominant in X-rays). 



averaged broad-band energetics, and provide a qualitative descrip- 
tion of the expected variability. But all agree that the origin and the 
details of the rapid variability patterns are likely to be complicated. 
The case of GX 339-4 may well be similar, and a full investigation 
of the parameter space of the various models is beyond the scope 
of this Letter. Nevertheless, several important conclusions can be 
drawn from our simultaneous multi-wavelength data. 



4.1 Reprocessing? 

The peak of the CCF time delay (150 ms) corresponds to a light- 
travel distance of 5000 Rg [= GM/c^] for A/>6 M©, too small 
for reprocessing on the companion star (GX 339-4 has a binary 
separation of ~ 25 light-seconds). An auto-correlation analysis 
of the individual light curves can also be used to constrain the 
emission processes. Fig. [3] shows the computed X-ray and optical 
auto-correlation functions (ACFs). Low count-rate Poisson noise 
dominating the X-ray ACF at zero lag has been corrected-for by 
subtracting white noise from an X-ray power spectrum, followed 
by an inverse Fourier transform. The final X-ray ACF is broader 
than the optical one, similar to the result found in the case of 
XTE Jl 11 8-1-480. This also argues against a reprocessing origin (on 
the outer parts of the accretion disk, say) for the rapidly- variable 
component, at least as described by a simple, linear transfer func- 
tion (cf. iKanbach et al. 2001) . 



4.2 Behaviour around flares and dips 

Fig. m shows that the behaviour of low- and high-energy X-ray 
(in this case, 2-5 keV and 5-20 keV) photon intensities is simi- 
lar, when selected around full band PCA peaks and dips, with no 
obvious lag. We note that CCFs of light curves extracted in these 
energy ranges, with respect to the optical, showed little difference 
to the full-band result of Fig.[T] The optical light curves follow the 
CCF shape, as already mentioned. No significant colour (g' — r') 
changes around the positions of X-ray flares/dips were detected. 

What about the source spectral behaviour? Using standard 
HEADAS tools, we extracted average X-ray spectra (and back- 
ground) within short time-bins, ~ ±50 — 150 ms, centred on 



Figure 4. Photon-index (F) contours for single power-law fits to the ex- 
tracted 3-30 keV flare and dip background-subtracted spectra, with absorp- 
tion fixed to Galactic. The total resultant time interval for spectral extraction 
is fa 30 s long for both flares and dips, and the net count rate is fa 40 and 5 
ct s~^ respectively. Contour levels correspond to 68.3, 95 and 99% for two 
interesting parameters, and include 1 per cent systematic errors. Results are 
shown for all three nights. The dotted diagonal line is the 1:1 line. 



X-ray PCA full-band flares and dips (HEXTE was not used for 
this analysis, as the source is always background-dominated above 
~ 20 keV). What we find is that the source hardens when it flares, 
and vice- versa. This is illustrated in Fig. [4] which shows the con- 
tours of independent fits to the extracted flares and dips spectra. A 
simple power-law (with absorption fixed to Galactic) was used to 
parametrize the spectral change between the two cases. The photon- 
index fitted to the flares spectrum is harder than that fitted to dips at 
99 per cent confidence on all nights. Simulations were used to con- 
firm that the results are robust to fluctuations in the background, 
which dominates above ^ 5 keV in the lower-flux dips spectrum. 

On the other hand, no significant changes were detected in X- 
ray spectra extracted at the positions of optical peaks as compared 
to optical dips - the spectral slopes were close to the slope for the 
average spectrum of the full dataset (F « 1.6). Extracting spectra 
150 ms before optical flares and dips (as suggested by the CCF 
delay) showed a small but clear difference in the intensities of the 
X-ray flare and dip spectra, but no obvious change in spectral slope 
within the errors. This seems consistent with the low CCF peak 
strength - i.e. every optical flare need not have been preceded by a 
locally-maximum X-ray flare, in spite of the average correlation. 



4.3 Implications 

The above analysis shows that the source is harder during X- 
ray flares, and vice-versa. In the context of a hot accretion disk 
corona model, this is consistent with increased Compton up- 
scattering when the source is (momentarily) brighter, and vice- 
versa. Fast Compton cooling timescales of tens of ms or smaller 
have been inferred in Galacti c black holes in the low/hard state 
(cf. iGuilbert & Fabianlll982h . for coronal electron temperatures 
kTe ^ 100 keV and seed-photon disk-blackbody X-ray lumi- 
nosities < 0.01I//LEdd- Very similar physical parameters are in- 
ferred for GX 339-4 du ring previous low/hard states ( Joi net et al.l 
120071 : iMiller et al.l 20061), and also contemporaneous with our ob- 
servations (cf. iTomsick et alj|2008l , who favour an inner disk ra- 
dius [rin] of ~ 10 Rg)- Thus, similarly-fast cooling timescales are 
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likely to hold for the case of GX 339-4, resulting in apparently- 
simultaneous flaring and hardening, which is what we observe. 

As for the optical emission, we have already excluded simple 
reprocessing models. Bremsstrahlung can also be ruled out, as this 
would require a corresponding X-ray flux higher than that observed 
by several orders of magnitude. The most likely remaining physical 
mechanism is then synchrotron emission. Cyclo-synchrotron mod- 
els in which several magnetized active regions with sizes of a few 
Schwarzschild radii contribute signi ficantly to the d ereddened op- 
tical flux have been investig ated by idi Matteo et al.l il99£, in ad- 
dition cf. iFabian et al.lll982h . The active regions are characterized 
by -B ~ few x 10^ G, coronal opti cal depths r ^ 0. 2 — 1 and 
fcTe ~ 150 - 200 keV. Fwardzihski & Zdziarskil feOOOh have also 
discussed important modifications to such models. 

Within the context of these models, what could the shape 
and the time delay of the CCF correspond do? The steep decline 
that follows the CCF peak (Fig. [TJ suggests the presence of some 
mechanism that cuts off the optical flares suddenly - like infall of 
synchrotron-emitting blobs into the BH, say. A potential difficulty 
of this scenario is that a free-fall time of 150 ms for a 6-10 M0 
BH corresponds to a physical radius of ~ 200 — 250 Rg, wh ich is 
larger than the value of rin inferred bv lTomsick et al.l 1 I2OO 8V Un- 
less the disk has receded further or an atypically large corona is 
present, this does not correspond to an obvious physically mean- 
ingful scale. 

An alternate hypothesis is that the rapid variability origi- 
nat es in non-thermal e mission within a relativistic outflow or a 
jet. iMalzac et alj ( 120041) invoked a magnetic energy reservoir that 
feeds both a jet component (dominating the optical), as well as 
an electron corona (dominating in X-rays) to explain the CCF of 
XTE Jl 118-1-480. Energy injection was modeled as shot flares. 
Feedback between the two components resulted in complex corre- 
lations, with the optical power being proportional to the differential 
of the X-ray flux. In such a model, the CCF shape and time lag are 
determined by the dissipation timescale of the process that injects 
energy into the jet. This should be largely independent of the exact 
injection mechanism, though the authors discuss the context of a 
magnetic energy reservoir. For our observations, a simple exponen- 
tial (oc e~'' ^) fit to the innermost part of the X-ray ACF (lags <0.5 
s) in Fig.[3]gives a dissipation timescale r = 0.2 ± 0.05 s (90% er- 
ror). This agrees with the observed optical delay, and suggests that 
such a differential correlation may apply to GX 339-4 as well. 

If X-ray flares above the accretion disk trigger the dissipation 
and large-scale re-ordering of the poloidal magnetic field thread- 
ing a jet, it is possible that the synchrotron (optical) emission will 
respond on timescales related to subsequent field build-up. Signifi- 
cant modulation of the poloidal field can occur on timescales orders 
of magnitude longer than the dynamica l time of the inner accretion 
disk regions where this field originates JLivio et alj2003l , see their 
Eq. 4). For rin ~ 10 -Rg, the Keplerian dynamical time is 6 ms. 
Our observed delay of the optical (i.e., jet) component {^^ 150 ms) 
is 25 times longer, which can easily be accommodated within the 
picture of magnetic modulation. Synchrotron emission by plasma 
accelerated along the jet during this period, followed by rapid ra- 
diative cooling, may thus explain the positive CCF peak and delay. 

As for the anti-correlation troughs: if the initial X-ray flares 
that triggered magnetic re-ordering are related to field reconnection 
events in a coexistent corona, the coronal magnetic energy density 
will be released on the X-ray flaring timescale. This will lead to 
a decrease of any ambient synchrotron emission that is occurring 
within the corona itself, resulting in an anti-correlation of optical 
with X-rays. X-ray flaring is coherent over timescales of ~ several 



seconds, as seen in the X-ray ACF (Fig. [3]l, and this is also the 
total length of the anti-correlation seen in the CCF plot of Fig.[T] In 
short, interaction between the jet and coronal components may be 
the key to understanding the complex correlation structure. 

A testable prediction of this model is that the appearance and 
strength of the positive CCF signal should be intimately related to 
the prominence of the jet. Our observations probed the beginning of 
the low/hard state. As the sources enters deeper into this state and 
the jet establishe s itself, its contribu tion to the overall energetics 
should grow (e.g. lFender et al.ll20041) . as should the related optical 
fractional variability rms. Support of this comes from the fact that 
we find an optical rms v ariability of ^ 15 p er cent over the full light 
curve timescales, while lMotch et al.l ( 119830 found an rms of 50 per 
cent when the source was optically brighter (this is w ell-matched 
to the estimate of ~ 50 per cent bv lCorbel & Fendenl2002 . based 
on broad-band photometry). Increased corona/jet coupling will re- 
sult in a higher peak strength of the optical/X-ray CCF. Finally, if 
a stronger poloidal field is required to establish a stronger jet com- 
ponent, this may result in longer timescales for breaking and estab- 
lishing this field component following reconnection flares. Chang- 
ing CCF delays could then be used to directly probe the evolution 
of characteristic accretion/ejection structures in X-ray binaries. 
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